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■ Abstract. The plethora of asymmetric planetary nebulae and the curi- 

ng \ ously high momenta of pre-planetary nebula outflows suggest that rota- 

^ 1 tional energy is extracted from the engines. Magneto-rotational outflows 

, driven by dynamos might be operating therein. I summarize scenarios 

' involving (a) an isolated AGB star and (b) a binary. The efficacy of (a) 

requires re-establishing differential rotation over the AGB star's lifetime, 
{Sj ' whereas (b) delivers the angular momentum at the end of the AGB phase 

when needed, and may involve turbulent disks. Both can produce fields 
that drive outflows of the required mechanical luminosity and momentum, 
though weak points and open questions require further study. 
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. 1. Observations Suggest a Role for Magneto- Rotational Outflows 

o _ 

^ <— I Powerful anisotropic outflows appear in a wide range of astrophysical sources 

JZLi' and magneto-rotationally driven models are often a favored explanation. Maybe 

asymmetric outflows (e.g. Balick & Frank 2003) in planetary nebulae (PPNe) 
£h ' and pre-planetary nebulae (PPNe) also involve dynamically important mag- 

c/2 . netic fields (Aller 1958). Magneto-rotational driving is particularly appealing 

for PPNe winds which are more powerful than PNe winds: The observation- 
ally inferred rates of mechanical momentum and energy deposition for PNe 
are (Kwok 2000) IT ~ lO 27 (M pn /O.O4M )(?;p ri /4Okm/s)(AtpA r /lO 4 yr)- 1 erg/cm, 
■ and E ~ 10 34 (II/10 27 erg/cm)(u pn /40km/s) erg/s, where M pn is the swept up 

mass in the PN, v pn is the speed of this mass, and the acceleration time Atpjy < 
age of the PN. But observations of PPNe (Bujarrabal et al. 2001) reveal val- 
ues up to 13 < 3 x lO 29 (M ppn /O.5M )(w ppn /2OOkm/s)(At ppn /lO 3 yr)- 1 erg/cm, 
and E < 10 37 (n/10 40 erg/cm)(u ppn /200km/s)erg/s where the acceleration time 
A.tpp n << age of the PPN. The PPNe wind momenta are often ~ 10 3 times that 
carried by radiation even when multiple scattering is considered. 

A large scale magnetic field can act as a drive belt that extracts rotational 
energy to drive a wind from the engine where the field lines are anchored (e.g. 
Spruit 1996). When the engine is a star, the extraction spins the star down. 
When the engine is an accretion disk, a fraction of the accretion luminosity gets 
redirected into the wind. In the Blandford & Payne (1982) type model, mass 
outflow is driven by the toroidal magnetic field pressure gradient along poloidal 
field lines, though in the rotator's co-rotating frame near the footpoints, the 
processes can be thought of as centrifugal launching along favorably inclined 
poloidal field lines. Consideration of magnetically mediated outflows raises nat- 
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ural questions in the PPNe/PNe context: Where do the required fields come 
from? How do these mechanisms work for single star or binary/disk systems? 
Are the outflow powers and momenta consistent with the observations? 

2. Why Would Dynamos be Needed to Produce the Fields? 

Convection driven and magneto-rotationally driven MHD turbulence (Balbus <fe 
Hawley 1998) are fully developed in in AGB stars and in accretion disks (the 
two relevant environments) respectively. Large scale magnetic fields can then 
experience turbulent diffusion so intense that the field is not flux-frozen on the 
relevant dynamical or advection time scales (Lubow et al. 1994; Blackman & 
Tan 2003 (BT)). This highlights the need for dynamo amplification. The fields 
also be large enough to avoid diffusing before escaping to the corona where they 
are needed. In coronae, the fields can further open up. Coronal loops which 
carry magneto-centrifugal winds should be at least as large as the radial scale 
of the anchoring rotator so that material can be accelerated to super-Alfvenic 
velocities before reaching the loop-tops and avoid sliding back down. 

3. Some Recent Developments in Large Scale Dynamo Theory 

Dynamos amplify or sustain magnetic fields by draining energy from random 
and shear motion. Nonhelical (direct) dynamos amplify fields at wave numbers 
k > kf whereas helical (inverse) dynamos can also amplify fields at k < kf, where 
kf is the smallest wavenumber at which the flow is turbulent. Here "helical" 
refers to whether the driving turbulence possesses a non-vanishing pseudoscalar 
helicity such as kinetic helicity (v • V x v), where v is the turbulent velocity, 
and the brackets imply averaging within a hemisphere. The nonhelical dynamo 
increases magnetic energy as the turbulence stretches magnetic field lines in a 
nearly random walk, but it is the helical dynamo (e.g. "a^fi," see Moffatt 1978) 
which produce the large scale fields needed for outflows. Helical turbulence can 
be supplied by convection in a star, or by the magneto-rotational instability in 
a stratified disk (Brandenburg et al. 1995; Brandenburg 1998; BT). 

While 3-D MHD simulations are essential for understanding dynamos, they 
are not always practical for modeling. For the latter, mean field theory (e.g. 
Moffatt 1978) is used. Here the field is broken into mean and fluctuating compo- 
nents. The helical dynamo growth coefficients become correlations of fluctuating 
quantities, and the mean field evolution is solved for. Traditional treatments do 
not dynamically determine the saturation strength of the large scale field; they 
either take the growth to be steady (=kinematic theory), or presume a form of 
quenching by the growing magnetic field without conserving magnetic helicity, 
a fundamental invariant of ideal MHD. A recent theory which alleviates these 
limitations by including the dynamical evolution of the velocity and magnetic 
helicity (Blackman & Field 2002) is illustrated below. 

Fig. 1 below (from Blackman &: Brandenburg 2003) shows the helical a^Tl 
dynamo in a rotating northern hemisphere. The kinematic theory is shown in 
(a) and (b) while (c), (d) and (e) show the dynamic theory, (a): To conserve 
angular momentum, rising (falling) blobs in a stratified medium twist in the 
opposite (same) sense to the underlying rotation. This implies a non-vanishing 
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(v • V x v), averaged in a given hemisphere. This "a^" effect implies that a large 
scale toroidal loop in the northern hemisphere incurs a right-handed writhe, and 
a radial field component as it rises, (b): Differential rotation (the "fT effect) at 
the base of the loop shears the radial component, amplifying the toroidal field. 
Dissipation or ballooning/opening of the top part of the loop allows for a net flux 
gain through the rectangle, (c): Same as (a) but now with the field represented 
as a ribbon. The right-handed writhe of the large scale loop is accompanied by a 
left-handed twist along the tube, incorporating magnetic helicity conservation, 
(d): Same as (b) but with field represented as ribbon/tube, (e): Top view of the 
combined twist and writhe The backreaction force which resists bending is the 
small scale twist. Diffusing the top part of the loops allows a net flux generation 
in the rectangles of (a)-(d), and helps alleviate the backreaction. 




The dynamic role of magnetic helicity is seen in simulations (Brandenburg 

2001) and modeled analytically (Field & Blackman 2002; Blackman & Field 

2002) . Although the small scale magnetic helicity growth drives the helical 
dynamo to a steady state, the large scale field grows large enough for most 
astrophysical applications. Application of these principles to disks and stars 
with boundaries is work in progress. In addition, buoyancy rather than kinetic 
helicity may be the initial driver of the helical a^ effect in disks (Brandenburg 
1998). This affects the initial sign of ay in each hemisphere, but magnetic helicity 
evolution remains a key to understanding the saturation. 



4. Single Star Scenario: AGB Dynamo and Magnetic Explosion 

This scenario (Blackman et al. 2001 (BFMTV)) begins with a main sequence 
(MS) star ~ 3M of rotational speed ~ 200 km/s. At the end the MS lifetime, 
the core collapses and the envelope expands, and angular momentum is assumed 
to be conserved on spherical mass shells. The rotation profile satisfies f2 oc 1/r 2 
which is then combined with a numerical AGB stellar model (SDK1401, Steve 
Kawalar, personal communication). In SDK1401, the convection zone extends 
from the outer stellar radius 5 x 10 12 cm down to the core radius r c = 9 x 10 10 cm. 
At the "interface" between the core and the envelope, Q ~ 2 x 10~ 5 cm/s. 

The combination of a strong differential rotation layer and overlying con- 
vection zone may lead to an "interface dynamo" (e.g. Markiel & Thomas 1999). 
Here the poloidal field is generated by a helical a^-effect layer in the convection 
zone. Diffusive transport (or turbulent pumping; Tobias et al. 1998) bleeds the 
poloidal field to the shear layer beneath the convection zone where it is further 
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amplified by linear stretching ("fT effect). BFMTV applied the code of Markiel 
k, Thomas (1999) to the AGB star, (The formula used was consistent with 
magnetic helicity conservation.) At r ~ r c the saturated field is B ~ 5 x 10 4 G. 

The scenario requires this dynamo to be steady for ~ 10 6 yr until the 
envelope of the AGB star is radiatively driven away in the late AGB stages. 
Eventually, the density falls such that the field can no longer be stored: Tur- 
bulent pumping can be gauged by assuming it to be proportional to a fixed 
convective energy density transfer rate pv 3 . As p decreases, v increases and 
the turbulent pressure, oc pv 2 , decreases with decreasing p, eventually drop- 
ping below B 2 /8it at the core. At this point, a "magnetic explosion" asso- 
ciated with toroidal pressure could ensue and drive the PPNe outflow. The 
magnetic explosion is driven by Poynting flux, and surface integration reveals 
a maximum mechanical luminosity (Blackman, Frank, Welch 2001 (BFW)) 



dius and angular speed respectively. This high L m lasts for a magnetic spin 
down time of the star (~ 100 yr, starting after the envelope is lost) and then de- 
cays exponentially. The asymptotic outflow speed should approach ~ SIota, 
where va is the Alfven radius, which can be less than or comparable to the escape 
speed v esc ~ 520 (r c /10 11 cm) 1 ^ 2 (M^/Mq) 1 ^ 2 km/s, depending on ta- Matt et 
al. (2003) (these proceedings) have simulated a magnetic explosion. The above 
results are consistent with the requirements of PPNe (see Sec. 1). (Imai et al. 
(2002) have observed a collimated AGB jet.) 

Unfortunately, the present interface dynamo models, while including a sat- 
uration of the ad effect, do not include saturation of the fi effect. In reality, 
amplification of the magnetic field at the shear layer drains shear energy, and 
transfers angular momentum from the core to the envelope, slowing down the 
core within < 100 yr, unless the differential rotation is re-seeded throughout 
the 10 6 yr lifetime of the AGB phase. This is important because matching the 
PPNe observations requires the field to be anchored in a rapidly rotating core. 

The so called "Lambda effect" (c.f. Riidiger 1989) might be able to re-seed 
the differential rotation. This mechanism involves the interplay between the 
overall rotation and asymmetric convection to produce a steady-state rotation 
profile. The helical interface dynamo has to be solved with a dynamically coupled 
Lambda effect. This awaits attention, but note that there is enough energy in 
the fusion driven turbulent convection to sustain the differential rotation: The 
turbulent energy density can be calculated from from the aforementioned AGB 
stellar model SDK1401 and the result is (J. Nordhaus personal comm.) ~ 5x 10 43 
erg. This can be compared with the magnetic field energy in the shear layer from 
BFMTV, which is E m ~ 10 41 erg, which in turn is about 1/8 the shear energy 
in the interface layer. Thus only ~ 2% of the energy from the convection needs 
to be steadily drained into the desired differential rotation. 



5. Binary Scenarios: Disk Dynamos and/or Primary Core Spin- up 

Binary scenarios do not suffer from the sustained shear problem just discussed 
because they deliver the angular momentum near the end of the AGB phase 
(e.g. Soker 1998), just when needed to drive the outflows. 




where r c and Cl c are the core ra- 
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Such angular momentum transfer could occur in the late AGB stage by- 
common envelope (CE) evolution (e.g. Iben & Livio 1993), as the secondary star 
penetrates the distended AGB envelope and transfers orbital angular momentum 
to it. The AGB star can shed 80% of its mass within years of the onset of 
CE, depending on the binary parameters. During the shedding, the companion 
spirals nearer to the core. While accretion of envelope material can occur onto 
the secondary (e.g. Morris 1987; Soker Sz Livio 1994; Mastrodemos Sz Morris 
1998), much larger (needed) accretion rates arise when an accretion disk forms 
around the primary via disruption of the secondary (Reyes- Ruiz &: Lopez 1999). 
Favored systems for this mode of accretion involve an evolved AGB star with 
mass 2.6 < M*/Mq < 3.6, a secondary with mass (< 0.1M@), and the initial 
binary separation > 200i? Q . Reyes- Ruiz & Lopez (1999) find that the disk 
accretion rate onto the residual primary core evolves in time (after an initial 
~ lyr viscous adjustment period) as = M ( i o (t/yr) _5 / 4 M yr _1 . Typical 
values of the initial accretion rate Mdo — lCT 3 M/yr. 

Once the disk forms, a helical disk dynamo (e.g. Pudritz 1981; Branden- 
burg et al. 1995; Reyes-Ruiz & Stepinksi 1995; Campbell 2000; BT) provides 
the mean field strengths from which one can determine the magnetic luminos- 
ity L m , available to drive an outflow. This can be estimated as the integrated 
Poynting flux flowing through the disk surface (or, the Poynting flux through the 
Alfven surface (BFW; Frank & Blackman 2003 (FB)). Tan & Blackman (2003) 

find L m ~ B p B^rf ~ ~ 10* (ff)^ (jfe) ( T ^)" (^f ? 

where a ss < 0.1 is the assumed disk viscosity parameter (Shakura & Sun- 
yaev 1979), V{ is the inner radius, and c s and J7j are the sound and rotation 
speeds there. The M& enters because p oc Md, and the mean surface poloidal 

and toroidal fields satisfy B p ~ p 1 / 2 a ss c s and Sa A calcula- 

tion of acceleration along a poloidal field line gives the asymptotic velocity 

Voo ~ SUr A ~ 2 1 / 2 ^r i ~ 730 (M*/M )5 (r i /10 11 cm)^ km/s (FB). For a rota- 
tionally supported disk, the asymptotic wind velocity is always near or greater 
than the escape speed. The mechanical luminosity and momentum delivered are 
consistent with the PPNe requirements of Sec. 1. 

Rather than form a disk as just described, an alternative is that the compan- 
ion could spin up both a slowly rotating core and envelope, leading to significant 
differential rotation between the two (assuming a turbulent viscosity in the en- 
velope). This could rejuvenate the interface dynamo of Sec. 4 during CE in the 
late AGB stage just when needed. The time scales of spin-orbit synchroniza- 
tion (Zahn 1977), envelope ejection, and viscous dissipation (Iben & Livio 1993) 
must be compared, but this alternative might mean that a wider range of binary 
systems could produce magneto-rotational outflows. 



6. Open Questions 

Asymmetric outflows in PPNe and PNe signature the transport of angular mo- 
mentum, and magnetically mediated outflows may produce the required power 
and momenta. Maser (Miranda et al. 2001; Vlemmings et al. 2003) and core 
X-ray (Kastner et al. 2003) observations loosely support this general paradigm. 
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However, key unresolved issues remain: (1) Is there a need for collimated 
outflows from both the AGB star and disk? Are both self-collimated or does 
the latter collimate the former? What might misaligned nested winds produce 
(BFW)? (2) Do all PPN/PNe produce collimated jets when the sources are 
young? (3) Large scale nonlinear dynamo theory is in its infancy for realistic 
boundary terms. How do helical dynamo fields open up and dynamically relax 
in the coronae? (4) How far is momentum carried as Poynting flux? (e.g. hydro- 
magnetic "fling" models vs. magnetically dominated "spring" models, or nested 
hybrid models? Ustyugova et al. 2000). (5) MHD jet simulations have been 
separate from MHD dynamo simulations. Simulations capturing both together 
have not been done, though using a dynamo produced field and then simulating 
the outflow provides encouraging results (von Rekowski et al. 2003). Non-steady 
calculations are important in this context. (6) The single AGB wind scenario will 
only work for PPNe if convection + rotation can steadily re-seed the differential 
rotation to ensure a large enough Poynting flux at the end of the AGB phase. 
If not, binaries may be absolutely required for asymmetric PNe and PPNe (e.g. 
Soker & Livio 1994; Soker 2001). Stellar evolution models which include mag- 
netic fields and rotation are needed. (7) To what extent are the strong outflows 
in PPNe correlated with the presence of binaries? Binaries with brown dwarfs 
or planets are not easily detected, but are favored for accretion around the AGB 
core. (8) How restrictive are the initial binary parameters for which common 
envelope evolution simulations (Demarco et al. 2003) would produce a suitably 
accreting disk to power PPNe? Is this relaxed by the scenario at the end of 
Sec. 5? (9) Consequences of internal dynamo produced fields for surface X-ray 
emission in AGB stars and disks need more study. (10) Could winds from the 
primary vs. secondary disks be distinguished by abundance differences (FB)? 
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